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Abstract 


2  Reductive  dechlorination  catalyzed  by  organohalidc-respiring  bacteria  is  often  considered 

3  for  remediation  of  non-aqueous  phase  liquid  (NAPL)  source  zones  due  to  cost  savings,  ease  of 

4  implementation,  regulatory  acceptance,  and  sustainability.  Despite  knowledge  of  the  key 

5  dechlorinators.  an  understanding  of  the  processes  and  factors  that  control  NAPL  dissolution  rales 

6  and  detoxification  (i.e.,  ethene  formation)  is  lacking.  A  recent  column  study  demonstrated  a  5- 

7  fold  cumulative  enhancement  in  tetrachloroethene  (PCE)  dissolution  and  ethene  formation  (Amos 

8  el  al.,  2009).  Spatial  and  temporal  monitoring  of  key  geochemical  and  microbial  (i.e.,  Geobacter 

9  lovleyi  and  Dehctlococcoides  mccartyi  strains)  parameters  in  the  column  generated  a  data  set 

10  used  herein  as  the  basis  for  refinement  and  testing  of  a  multiphase,  compositional  transport 

1 1  model.  The  refined  model  is  capable  of  simulating  the  reactive  transport  of  multiple  chemical 

12  constituents  produced  and  consumed  by  organohalide-respiring  bacteria  and  accounts  for 

13  substrate  limitations  and  competitive  inhibition.  Parameter  estimation  techniques  were  used  to 

14  optimize  the  values  of  sensitive  microbial  kinetic  parameters,  including  maximum  utilization 

1 5  rates,  biomass  yield  coefficients,  and  endogenous  decay  rates.  Comparison  and  calibration  of 

!  6  model  simulations  with  the  experimental  data  demonstrate  that  the  model  is  able  to  accurately 

1 7  reproduce  measured  effluent  concentrations,  while  delineating  trends  in  deehlorinator  growth 

18  and  reductive  dechlorination  kinetics  along  the  column.  Sensitivity  analyses  performed  on  the 

19  optimized  model  parameters  indicate  that  the  rates  of  PCE  and  cis- 1 ,2-dichloroethene  (c/s-DCE) 

20  transformation  and  Dehalococcoides  growth  will  govern  bioenhanced  dissolution,  as  long  as 

21  electron  donor  (i.e.,  hydrogen  flux)  is  not  limiting.  Dissolution  enhancements  were  shown  to  be 

22  independent  of  c/.v-DCE  accumulation;  however,  accumulation  of  cA-DCE,  as  well  as  column 

23  length  and  flow'  rate  (i.e.,  column  residence  time),  strongly  influenced  the  extent  of  reductive 
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dechlorination.  When  cw-DCE  inhibition  was  neglected,  the  model  over-predicted  ethene 

2  production  ten-fold,  while  reductions  in  residence  time  (i.e.,  a  two- fold  decrease  in  column 

3  length  or  two-fold  increase  in  flow  rate)  resulted  in  a  more  than  70%  decline  in  ethene 

4  production.  These  results  suggest  that  spatial  and  temporal  variations  in  microbial  community 

5  composition  and  activity  must  be  understood  to  model,  predict,  and  manage  bioenhanced  NAPL 

6  dissolution. 

7 

8  Keywords :  PCE;  NAPL;  enhanced  dissolution;  modeling;  model  verification;  reductive 

9  dechlorination 

10 

1 1  1.  Introduction 

12  Chlorinated  solvents  (e.g.,  tetrachloroethene  [PCE],  trichloroethene  [TCE]),  are  widely 

13  used  in  industrial  processes,  and  at  military  installations  and  government  facilities,  and  have 

14  often  been  released  as  dense  non-aqueous  phase  liquids  (DNAPLs).  DNAPL  dissolution, 

15  whether  from  pools  retained  above  or  within  low-permeability  strata  or  from  residual  saturation 

16  ganglia  trapped  in  the  porous  soil  or  rock  matrix,  creates  a  persistent  source  of  groundwater 

1 7  contamination.  In  the  past  two  decades,  numerous  remediation  technologies  have  been  developed 

18  to  treat  DNAPL  source  zones;  however,  the  ability  of  a  single  technology  to  completely  remove 

19  or  destroy  all  DNAPL  mass  and  reduce  dissolved-phase  contaminant  concentrations  below 

20  drinking  writer  standards  is  limited.  Among  potential  in  situ  remediation  technologies,  microbial 

21  reductive  dechlorination  has  emerged  as  an  attractive  plume  remedy,  and  an  approach  to  enhance 

22  rate-limited  contaminant  mass  transfer  at  sites  where  contaminants  arc  sequestered  in  low 

23  permeability  zones  (e.g.,  Scheutz  el  al.,  2010)  or  present  in  DNAPL  source  zones  (Da  Silva  et  a!.. 
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1  2006;  Sleep  et  al,,  2006;  Schaefer  et  ah,  2010).  In  these  applications  microbial  reductive 

2  dechlorination  may  also  serve  as  a  source  zone  polishing  step  to  control  residual  contaminant 

3  concentrations  following  aggressive  physicochemical  treatment  (Mravik  et  al.,  2003;  Ramsburg 

4  el  al.,  2004;  Christ  et  al.,  2005).  During  DNAPL  source  zone  bioremediation,  microbial  activity 

5  acts  to  lower  dissolved-phase  contaminant  concentrations,  thereby  increasing  the  driving  force 

6  for  contaminant  dissolution  from  the  nonaqueous  to  the  aqueous  phase,  a  process  commonly 

7  referred  to  as  bioenhanced  dissolution  (Yang  and  McCarty.  2000;  Cope  and  Hughes,  2001 ;  Yang 

8  and  McCarty,  2002;  Adamson  et  al.,  2003:  Sleep  et  af,  2006;  Glover  et  al.,  2007;  Amos  et  al., 

9  2008).  While  dissolution  enhancements  can  lead  to  short  term  increases  in  dissolved  phase 

10  transformation  products  (i.e.,  cA-DCE  concentrations),  the  overall  source  longevity  and 

1 1  associated  cleanup  times  decrease,  potentially  resulting  in  reduced  long-term  risk. 

12  Results  obtained  from  batch,  column,  and  aquifer  ceil  experiments  indicate  that  the  extent 

13  of  DNAPL  dissolution  enhancement  due  to  bioactivity  can  vary  widely,  from  a  negligible  effect 

14  to  a  more  than  20-fold  increase  (e.g.  Carr  et  af,  2000;  Yang  and  McCarty,  2000;  Cope  and 

15  Hughes,  2001;  Sleep  et  af,  2006;  Amos  et  af,  2008,  2009;  Schaefer  et  al.,  2010;  Philips  et  al., 

16  2011).  This  variability  in  dissolution  enhancement  has  been  attributed  to  several  factors, 

17  including  differing  microbial  populations,  electron  donor  limitations,  competitive  and  threshold 

1 8  inhibitions,  and  competitor  populations  (e.g.,  methanogens),  as  well  as  changing  flow  conditions 

19  due  to  gas  formation,  microbial  growth,  and  experimental  design.  While  bench-scale 

20  experiments  provide  excellent  demonstrations  of  the  conditions  leading  to  bioenhanced  DNAPL 

21  dissolution  (e.g..  Amos  et  al.,  2009;  Schaefer  et  af,  2010),  their  use  as  a  design  tool  for  treatment 

22  applications  in  the  field  is  limited.  Numerical  simulators  provide  a  complementary  tool  that 

23  allows  for  more  efficient  assessment  of  individual  parameters  and  inter-related  processes 
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impacting  microbial  activity,  guiding  future  efforts  to  understand  and  optimize  in  silu 
bioenhanced  dissolution. 

Prediction  of  bioenhanced  DNAPL  dissolution  requires  a  model  that  quantifies  inter¬ 
relationships  between  dechlorination  and  mass  transfer  processes.  A  number  of  models  capable 
of  simulating  organohalide  respiration  kinclics  have  been  presented  in  the  literature  that 
incorporate  electron  donor  flux  (i.e.,  fermentation)  and  competition  (i.e.,  methanogenesis)  {e.g. 
Bagley,  1998;  Fennel!  and  Gossett,  1998;  Lee  ct  al.,  2004;  Yu  et  al.,  2004,  2005,  Becker,  2006; 
Huang  and  Becker,  2009;  Becker  and  Seagren,  2009;  Haest  et  al.,  2010;  Kouznetsova  et  al.,  2010; 
Sabalowsky  and  Semprini,  2010;  Huang  and  Becker,  2011).  These  models,  however,  are  limited 
to  the  simulation  of  dechlorination  in  batch  reactors  or  single-phase,  aqueous  systems.  While 
several  models  that  include  both  dechlorination  and  dissolution  have  been  presented  (Chu  et  al., 
2003;  Widdowson,  2004;  Hammond  et  al..  2005:  Chambon  et  al.,  2010),  those  models  have  only 
been  applied  to  simplified  systems  (e.g.,  no  dissolution)  and  thus,  have  not  considered  the 
potential  influence  of  a  diminishing  DNAPL  mass  on  the  temporal  and  spatial  evolution  of  the 
mobile  aqueous  phase  flow’  and  microbial  activity.  Although  a  number  of  multi-phase,  multi- 
component  simulators  have  been  adapted  to  model  organohalide  respiration  in  DNAPL  source 
zones  (Delshad  et  al.,  1996;  Abriola  et  al.,  1997;  Rathfelder  et  ah,  2000;  Willis  and  Shoemaker, 
2000;  Battisstelli,  2004;  Clement  et  ah,  2004),  very  few  studies  (Christ  and  Abriola.  2006; 
Becker  and  Seagren,  2009)  have  explored  dissolution  enhancements  due  to  microbial  activity. 
Christ  and  Abriola  (2006)  adapted  a  multiphase  compositional  simulator  (e.g,,  MISER  - 
Michigan  Subsurface  Environmental  Remediation  Simulator)  to  incorporate  dynamic  interphase 
mass  transfer,  non-linear  dechlorination  kinetics,  microbial  inhibition,  and  competition  in  a 
framework  that  is  capable  of  simulating  two-dimensional  non-uniform  source  zones.  Their 
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1  comparisons  to  available  data  indicated  that  the  model  was  capable  of  capturing  effluent 

2  concentrations  consistent  with  bioenhanced  dissolution.  However,  these  data  were  limited  to  the 

3  PCE-to-ds-DCE  dechlorination  step,  and  did  not  provide  the  information  necessary  to  fully 

4  calibrate  and  validate  the  model.  The  resulting  simulations  were  thus  hypothetical  in  nature  and 

5  provided  limited  insight  into  the  fundamental  factors  controlling  bioenhanced  dissolution 

6  efficacy  within  a  source  zone,  including  conditions  leading  to  incomplete  or  stalled  reductive 

7  dechlorination,  a  frequent  limitation  of  bioremediation  during  field  application. 

8  A  continuous-flow  column  experiment  (Amos  et  ah,  2009),  instrumented  to  quantify 

9  aqueous  phase  constituent  concentrations  (e.g.,  organic  acids.  PCE,  TCE,  czs-DCE,  vinyl 

10  chloride  [VC],  ethene)  and  the  distributions  of  key  dechlorinating  populations  (i.e.,  Geobacter 

1 1  lovleyi  strain  SZ  [GEOSZ]  and  Dehalococcoides  mccartyi  [DHC])  along  the  length  of  the 

12  column,  provides  a  novel,  high  resolution  data  set  to  elucidate  factors  controlling  bioenhanced 

13  dissolution  and  reductive  dechlorination  within  DNAPL  source  zones  and  associated  down- 

14  gradient  plume  regions.  In  this  work,  the  experimental  results  of  Amos  et  al.  (2009)  are 

15  compared  to  model  simulations  to:  (i)  examine  the  ability  of  a  state-of-the-art  mathematical 

16  model  to  simulate  spatial  and  temporal  changes  in  bioenhanced  dissolution  and  reductive 

17  dechlorination;  (ii)  quantify  microbial  transformation,  utilization,  and  growth  rates  using 

1 8  parameter  estimation  techniques:  and  (iii)  identity  and  assess  system  parameters  most  influential 

19  in  controlling  bioenhanced  dissolution  and  reductive  dechlorination  within  a  DNAPL  source 

20  zone.  Knowledge  gained  from  the  analysis  of  the  relative  importance  of  system  parameters  to 

21  each  step  of  the  dechlorination  process  provides  additional  insight  into  dechlorination  and 

22  bioenhanced  dissolution  mechanisms  and  may  serve  to  direct  future  research,  and  guide  effective 

23  implementation  of  bioenhanced  dissolution  strategics  at  the  field  scale. 
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]  2.  Methodology 

2  2.1  Column  experiment 

3  Bioaclive  column  data  were  primarily  obtained  from  Amos  et  al.  (2009),  in  which  a  4.8 

4  cm  (inside  diameter)  x  60  cm  (length)  borosilicate  glass  column  was  constructed  with  1 1  glass 

5  sampling  ports  located  on  alternating  sides  of  the  column  at  5  cm  intervals  to  allow  for 

6  monitoring  of  aqueous  phase  constituents  and  microbial  populations  in  space  and  time  (Figure  1 ). 

7  The  column  was  packed  with  sterile  Federal  Fine  Ottawa  sand  under  anoxic,  water-saturated 

8  conditions.  The  column  was  uniformly  inoculated  with  a  non-methanogenic,  PCE-to-ethene 

9  dechlorinating  microbial  consortium  (BDI-SZ)  that  contained  DHC,  a  Dehalobacter  sp.  (DUB), 

10  and  GEOSZ.  The  porosity  of  the  packed  column  was  approximately  0.37  (Table  1),  yielding  a 

1 1  total  pore  volume  of  407  mL.  During  the  first  stage  (days  1  -  4)  of  the  experiment,  a  reduced 

12  mineral  salts  medium  without  growth  substrates  was  flushed  through  the  column  in  an  up  flow 

13  mode  at  a  rate  of  0.25  mL/min  for  3.3  pore  volumes.  Following  this  microbial  elution  phase,  a 

14  10-cm  long  NAPL  source  zone  was  established  by  injecting  an  anoxic  NAPE  containing 

15  0.25:0.75  (mol: mol)  PCE  in  hexadecane  into  the  first  10-15  cm  of  the  column  (upflow  mode), 

1 6  after  which  the  flow  direction  was  reversed  and  a  reduced  medium  was  introduced  from  the  top 

17  of  the  column  (downflow  mode)  to  displace  mobile  NAPL  from  the  pore  space.  The  resulting 

18  source  zone  area  contained  approximately  14  mL  of  PCE-eontaining  NAPL  that  was  uniformly 

19  distributed  over  the  first  10  cm  of  the  column  as  entrapped  droplets  and  ganglia  (average  initial 

20  saturation  =  0.21),  which  transitioned  (10-15  cm)  to  a  NAPL-free,  down-gradient  plume  region 

2 1  over  the  remainder  of  the  column  (15-60  cm). 

22  Following  establishment  of  the  NAPL  source  zone,  a  reduced  medium  amended  with  20 

23  mM  lactate  was  delivered  in  an  upflow  mode  at  a  flow  rate  of  0.25  mL/min,  corresponding  to  a 
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1  residence  time  of  1.1  days.  Following  the  injection  of  16.8  pore  volumes  (19  days),  a  22  hour 

2  flow  interruption  occurred  (caused  by  an  equipment  malfunction).  The  flow  then  resumed  at  the 

3  same  flow  rate,  and  after  21 .6  pore  volumes  (24  days),  it  was  reduced  to  0.1  mL/rnin  to  increase 

4  the  solution  residence  time  within  the  column.  Effluent  samples  were  collected  every  1  to  3  days 

5  and  side-port  samples  were  collected  periodically  (after  9,  14,  20,  28  and  30  pore  volumes)  to 

6  monitor  for  chlorinated  ethenes,  ethene,  pH,  organic  acids,  and  dechlorinating  bacteria.  After 

7  PCE  was  completely  depleted  from  the  mixed-NAPL  (ca.  day  55),  a  1.1  pore  volume  pulse  (450 

8  mL)  of  reduced  medium  amended  with  VC  was  injected  into  the  column.  This  final  stage  of  the 

9  experiment  was  used  to  determine  if  the  BDI-SZ  consortium  was  capable  of  complete 

10  dechlorination  of  VC  to  ethene  following  depletion  of  polychlorinated  ethenes. 

1 1  During  the  first  4  days  of  NAPL  dissolution,  effluent  PCE  concentrations  were  near 

12  equilibrium  (Figure  2),  which  is  consistent  with  the  delay  In  the  onset  of  lactate  fermentation, 

13  suggesting  that  PCE3  dechlorination  was  limited  by  the  absence  of  direct  electron  donors  (i.e., 

14  hydrogen  and  acetate).  Between  days  4  and  1 2.  the  total  chlorinated  ethenes  concentration  (PCE. 

15  cfv-DCE  and  VC)  approached  the  system  equilibrium  value,  indicating  that  while  dechlorination 

16  was  occurring,  enhancements  in  dissolution  due  to  bioactivity  were  negligible.  Amos  et  al. 

1 7  (2009)  hypothesized  that  this  phase  of  the  experiment  represented  an  acclimatization  period,  in 

18  which  the  microorganisms  reached  sufficient  numbers  in  the  NAPL-contaminated  source  zone 

19  region  to  facilitate  enhanced  dissolution.  After  day  12,  a  rapid  increase  in  efv-DCE  formation 

20  was  observed,  consistent  with  bioenhanced  dissolution.  A  spike  in  VC  was  observed  after  the 

21  flow  interruption  on  day  19,  but  VC  concentrations  returned  to  previous  levels  when  flow 

22  resumed.  This  observation  suggests  that  the  column  residence  time  was  insufficient  for  the  cis- 

23  DCE  to  VC  dechlorination  step.  After  the  flow  rate  was  reduced  from  0.25  to  0.1  rnL/min  on 
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1  day  24,  VC  concentrations  continued  to  increase  until  day  45,  when  the  PCE  in  the  column  was 

2  completely  depleted  and  ethene  production  began.  The  observation  that  ethene  formation  was 

3  minimal  until  after  day  45  suggests  that  VC  to  ethene  dechlorination  was  inhibited  by  high  cis- 

4  DCE  concentrations  (e.g.,  Yu  ct  ah,  2005).  These  results,  when  combined  wath  side-port  sample 

5  data,  provide  a  comprehensive  data  set  to  support:  (i)  conceptual  and  mathematical  model 

6  development  and  validation  and  (ii)  the  subsequent  investigation  of  the  influence  of  residence 

7  time  and  c7.s-DCE  formation  on  complete  detoxification  of  PCE  to  ethene  in  a  NAPE  source  zone 

8  and  associated  plume  region. 

9  2.2  Model  development 

10  Numerical  simulations  were  performed  using  a  modified  version  of  MISER  (Christ  and 

1 1  Abriola.  2006),  a  two-dimensional,  finite  element  immiscible  multiphase,  multi-component 

12  transport  simulator  originally  developed  for  soil  vapor  extraction  and  bioventing  (Abriola  et  al., 

13  1997:  Rathfelder  et  al..  2000).  Christ  and  Abriola  (2006)  modified  the  original  biodegradation 

14  module  in  MISER  to  incorporate  metabolic  reductive  dechlorination  of  multiple  contaminants  by 

1 5  multiple  microbial  populations  in  a  DNAPL  source  zone.  The  detailed  mathematical  formulation, 

16  including  dechlorination  kinetics,  and  its  numerical  implementation  have  been  presented 

17  previously  (Abriola  et  al..  1997;  Christ  and  Abriola,  2006).  The  main  features  of  the  modified 

1 8  MISER  model  are  described  below. 

19  2.2 A.  Governing  equations 

20  The  NAPI.  phase  is  assumed  to  be  discontinuous  and  immobilized  by  capillary  forces, 

21  while  the  aqueous  phase  is  assumed  to  be  mobile  and  influenced  by  changes  in  relative 

22  permeability  as  the  NAPE  dissolves.  Phase  mass  balance  equations  are  given  as: 
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3  where.  qu  =$SWVW  = - —  (VPH.  -pwg)  is  the  aqueous  phase  flux  computed  by  the  modified 

Mil 

4  form  of  Darcy's  law,  <fi  is  the  matrix  porosity  (-),  p,,  and  p0  are  the  aqueous  and  NAPL  phase 

5  density  (M/L3),  Sw  and  S0  arc  the  aqueous  and  NAPL  phase  saturation  (-),  and  ELwa ,  BLlv , QK  are 

6  the  interphase  mass  exchange,  bioreaction,  and  source/sink  terms  (M/L’-t),  respectively.  Eq.  (la) 

7  relates  the  change  in  aqueous  phase  mass  to  the  aqueous  phase  advective  flux  (first  term),  mass 

8  transfer  to  the  organic  phase  (second  term),  biological  reactions  (third  term),  and  source/sinks 

9  (fourth  term).  Mass  transferred  from  the  aqueous  phase  accumulates  in  the  organic  phase  as 

1 0  shown  in  equation  1  b. 

1 1  Within  the  aqueous  phase,  the  mass  balance  for  each  chemical  component  c  is  expressed: 

12  *^SJZ)+0-S.(v&-D‘.'7C:)-BU+K  (2) 

13  where  C‘.  is  the  concentration  of  component  c  in  the  aqueous  phase  (M/I.3),  v„  is  the  pore  water 

14  velocity  (L/t),  D‘,  is  the  aqueous  phase  hydrodynamic  dispersion  coefficient  for 

1 5  component  c  (L  /t),  and  all  other  parameters  are  as  given  previously.  Note  that  since  no  mass 

16  transport  occurs  within  the  NAPL  phase,  there  is  no  corresponding  equation  for  immobile  NAPL. 
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1  2.2.2  PCE  Dissolution 


2  PCE  mass  transfer  from  entrapped  NAPL  to  the  aqueous  phase  (i.e.,  dissolution)  is 

3  modeled  using  a  linear  driving  force  expression  to  approximate  the  concentration  gradient 

4  between  the  surface  of  the  NAPL  and  the  bulk  aqueous  phase: 

5  £”=e"(cr--c*)  o) 

6  where  C'JJ'"''  is  the  aqueous  concentration  in  equilibrium  with  the  mixed  NAPL,  as  determined 

7  by  Raoulfs  law  and  the  assumption  of  ideal  fluid  behavior  (i.e.,  activity  coefficients  equal  one) 

8  (Schwarzenbach  et  ah,  2003),  and  C[ih  is  the  bulk  aqueous  phase  PCE  concentration.  The 

9  lumped  mass  transfer  coefficient.  k’JJ''  ,  is  the  product  of  the  mass  transfer  coefficient  and 

10  interfacial  area  per  unit  volume  of  porous  media.  This  lumped  mass  transfer  coefficient  is 

1 1  typically  modeled  using  a  Sherwood  number  correlation  (e,g..  Miller  et  ah.  1990;  Powers  et  ah. 

12  1992,  Imhoff  et  ah,  1994).  Given  the  conditions  of  the  column  experiment,  the  modified 

13  Sherwood  number  correlation  for  transient  dissolution  developed  by  Powers  et  ah  (1994)  was 

14  judged  most  appropriate: 

r  \  0 .673  /  ^  xdr 

15  Sh  =  4. 1 3 Rea™  ^ -  U°m  •%  (4) 

UJ  UJ 

1 6  where  the  Reynolds  number  is  defined  as  Re  =  qwpwd50  /  //„, ,  di{)  / dM represents  the  median  grain 

17  size  normalized  by  a  “medium”  sand  grain  diameter  (d,u  )  of  0.05  cm.  as  defined  by  the  U.S. 

18  Department  of  Agriculture  (Driscoll,  1986),  Ut  =db0/dw  is  the  uniformity  index,  and  the  fitting 

19  parameter  a  =  0.518  +  0.1 14(t/5H  / dM )  +  0.  \Ut .  I  lere,  dj  denotes  the  grain  size  that  is  i  %  of  finer 


1 1 


1  particles  by  weight,  and  5,,  and  S0  are  the  initial  and  transient  NAFL  saturations,  respectively. 


2  Rearranging  the  modified  Sherwood  number,  Sh  =  /  Dt  (Miller  el  ah,  1990)  gives: 


kr<li-4]3^^Reom  —  1  uom 


4  Note  that  k'JJ  is  a  function  of  the  Reynolds  number  and  hence,  of  the  aqueous  phase  velocity. 


5  During  flow  interruption,  the  velocity  is  zero  and  Eq.  (5)  is  no  longer  applicable.  In  this  case,  the 

6  lumped  mass  transfer  coefficient  is  assumed  to  be  a  constant,  estimated  by  fitting  measured  data. 


7  2.2.3  Kinetic  model  for  reductive  dechlorination 

8  The  conceptual  model  presented  in  Christ  and  Abriola  (2006),  which  incorporates  PCE  to 

9  c/.v-DCE  (via  TCE  as  intermediate)  followed  by  cw-DCE  to  ethene  (via  VC  as  intermediate) 

10  reductive  dechlorination  reactions,  was  adapted  for  the  BDI-SZ  column  experimental  conditions. 

1 1  The  PC  E -to-ethene -deehlorinating  consortium  BDI-SZ  contains  populations  (FER)  that  convert 

12  lactate  to  acetate  and  propionate  and  increase  hydrogen  flux,  thus  providing  direct  electron 

13  donors  (i.e.,  acetate  and  hydrogen)  for  the  organohalide-respiring  populations  responsible  for 

14  PCE-to-ethene  reductive  dechlorination.  Methanogcns  are  absent  from  the  BDI-SZ  consortium 

15  ({Amos  et  al .,  2007b;  Amos  et  al .,  2008;  Amos  el  al.,  2009;  Ritalahti  el  al.,  2006)  and  thus, 

16  competitor  populations  were  not  considered  in  this  application. 

1 7  Chlorinated  solvent  degradation  was  assumed  to  follow  Monod  kinetics  (Fennell  and 

1 8  Gossett,  1 998;  Lee  et  al.,  2004;  Christ  and  Abriola,  2006); 
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where  model  parameters  are  defined  in  Table  2.  Note  that  an  inhibition  factor  f  is  included  for 
the  2lltl  dechlorinator  population,  Xdhc<  to  account  for  a  decrease  in  dechlorinating  activity  due  to 
the  presence  of  PCE  and  competition  among  dechlorinating  species  for  electron  acceptors  (Christ 
and  Abriola,  2006): 


r  = 


1,  c  =  DCE, 


!  + 
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jxi; 


c  =VC 


(7) 


7  Here,  k\  is  the  inhibition  coefficient  due  to  the  presence  of  c  =  PCE  and  cw-DCE.  As  shown  in 

8  Eq.  (7),  PCE  and  cA-DCE  inhibit  degradation  of  VC,  and  hence  the  production  of  ethene. 

9  Equations  (6)  and  (7)  can  be  combined  to  form  the  bioreaction  terms 7?‘(  in  cq.  (la)  and  eq.  (2) 

10  according  to:  ,  B1™  =-R’™  ,  Blvli  =  R'™-R™  ,  B1™  =  R?'«  - R!fH  , 

11  B[{  =R!nKh  -  R'J  ,  B!n 1 H  =  R'J  ,  and  B^1  =  F!lM  R'f  -  FR"r  ,  where  F"  is  the  stoichiometric 

t 

12  production  or  use  of  hydrogen. 

13  Microbial  growth  during  each  dechlorination  step  (  dX  l  dt  )  was  modeled  using 

14  microorganism  yield  (T’)  resulting  from  organohalide  respiration,  and  the  first-order  endogenous 

1 5  decay  rate  (kj)  {Yu  et  al.,  2005): 
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(8b) 


+  Ym:R™'H -kf°Xr,m 

dt 

2  =  Ynn  R'Jn  +  YVCR«'  -kfcXim,  (8c) 

dt 

3  where  Xj  is  the  attached  active  microbial  cell  concentration  and  kj  is  the  endogenous  decay  rate 

4  of  population  j.  Although  some  authors  have  proposed  that  kj  can  increase  due  to  inhibitory 

5  effects  (Sabalowsky  and  Semprini,  2010),  kj  was  assumed  constant  in  this  work  based  on  recent 

6  data  suggesting  inhibition  is  related  to  reduced  dechlorination  activity,  rather  than  enhanced  cell 

7  decay  (Philips  et  a!.,  2013).  The  D1IB  population  was  ignored  in  the  model  based  upon  the 

8  experimental  observation  that  DHB  did  not  colonize  the  column  and  hence  did  not  contribute  to 

9  the  dechlorination  kinetics  (Amos  et  ah,  2009). 

10  The  set  of  equations  (1)  through  (8)  was  implemented  in  MISER  using  a  Galerkin  Unite 

1 1  element  discretization  method  (Zienkiewicz  and  Taylor.  1991).  The  flow  and  transport  equations 

12  in  MISER,  along  with  the  constitutive  relationships,  have  been  verified  by  mass  balance 

13  calculations  and  with  analytical  or  other  numerical  solutions  (Abriola  ct  ah,  1997,  Rathfeldcr  et 

14  ah,  2000;  Christ  and  Abriola,  2006).  In  the  following  sections,  the  experimental  results  of  Amos 

15  et  al.  (2009)  are  used  for  model  calibration  and  to  evaluate  the  ability  of  the  model  to  reproduce 

16  observed  spatial  and  temporal  trends  related  to  the  reductive  dechlorination  of  chlorinated 

1 7  elhcnes  and  ethene  formation. 

18 

19  2.2.4  Parameter  Estimation 

20  Application  of  any  complex  multi-phase,  multi-constituent  mode! 

21  specification  of  a  large  number  of  system  parameters  and  initial  conditions. 


requires  the 
This  section 
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1  describes  the  approach  used  to  estimate  appropriate  model  parameters  and  any  associated 

2  changes  necessary  to  implement  these  in  the  simulator. 

3  2.2.4.  J  Initial  hydrogeologic  and  biomass  parameters 

4  As  described  in  Section  2.1  and  shown  in  Figure  2,  microbial  dechlorination  was  not 

5  observed  during  the  first  4  days  of  the  experiment,  and  dechlorination  that  did  occur  from  days  4 

6  to  12  did  not  enhance  PCE  dissolution.  One  possible  explanation  for  this  lag  is  the  organisms’ 

7  acclimation  period  (Amos  el  al.,  2009).  To  capture  this  effect  in  the  model,  the  maximum 

8  utilization  rates  for  dechlorination  activity  were  set  to  zero  for  day  0  to  day  4  and  the  initial  DHC 

9  biomass  concentrations  measured  in  the  two  source  zone  ports  (jc  <  0. 3  m)  were  replaced  with  the 

10  biomass  amounts  measured  in  those  ports  at  day  16  of  the  experiment;  the  first  available  side 

11  port  sample  data  (see  Figure  3  for  biomass  profiles).  While  some  growth  may  have  occurred 

12  during  the  first  4  days  of  the  experiment,  it  was  assumed  to  be  negligible,  consistent  with  this 

13  lack  of  dechlorination  activity.  Furthermore,  initial  numerical  experiments  suggested  that  the 

14  initial  biomass  concentrations  (day  0)  measured  by  Amos  el  al.  (2009)  were  likely  over-estimates 

15  of  viable  biomass,  which  was  attributed  to  the  introduction  of  PCE-containing  NAPL  in  the 

16  source  zone  region.  Initial  hydrogeologic  parameters  were  obtained  from  direct  experimental 

17  observations  (Table  1).  The  lumped  mass  transfer  coefficient  during  flow'  interruption  was 

18  estimated  to  be  i  .Ox  [  0“4  s'1  using  a  trial  and  error  fitting  process. 

19  2.2. 4. 2.  Biokinetic  parameter  estimates 

20  Microbial  kinetic  parameters  were  not  directly  available  from  prior  experimental  work, 

21  and  therefore,  were  estimated  using  inverse  modeling.  Estimates  were  needed  for  a  total  of  20 

22  biological  parameters:  the  maximum  utilization  rates  (k^m  ),  half-saturation  constants  (  K's)  and 
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1  yield  coefficients  ( Yc )  for  lactate  fermentation  and  reductive  dechlorination  of  PCE,  TCE,  cis- 

2  DCE,  and  VC  {15  parameters);  the  endogenous  decay  rate  ( ktl )  for  three  relevant  microbial 

3  groups  (FER,  D!  1C.  GEOSZ);  and  the  inhibition  coefficients  (ft,1  )  for  PCE  and  c/s-DCE  (2 

4  parameters)  {not  TCE  since  TCE  did  not  accumulate).  Initial  values  of  these  microbial  kinetic 

5  parameters  were  derived  from  the  literature  (Table  2).  Initial  model  calibration  efforts  revealed 

6  that  the  calibration  process  (inverse  model  convergence)  proceeded  more  rapidly  when  the  initial 

7  maximum  utilization  rates  obtained  from  the  literature  were  decreased  by  factors  of  four  to  20 

8  (Table  2).  The  need  to  adjust  biokinetic  parameters  was  expected  given  the  variability  of 

9  reported  dechlorination  rates  for  different  microbial  consortia  {Table  2).  Note  that  al!  parameter 

10  adjustments  fell  within  the  range  of  values  reported  in  the  literature. 

1 1  The  large  number  of  biokinetic  parameters,  coupled  with  the  correlation  between 

1 2  parameters  and  insensitivity  of  simulations  to  some  parameters,  resulted  in  an  ill-defined  matrix 

13  for  inverse  modeling  as  discussed  below.  Thus,  the  biokinetic  calibration  problem  was 

14  simplified  by  identifying  the  eight  most  sensitive  parameters:  maximum  utilization  rates  for  PCE, 

15  TCE,  cfv-DCE,  and  VC  (four  parameters),  biomass  yield  from  PCE  and  cis- DCE  degradation 

16  (two  parameters),  and  endogenous  decay  rates  of  GEO  and  DHC  (two  parameters).  Inverse 

17  optimization  efforts  were  focused  on  these  parameters,  which  were  considered  adjustable 

18  parameters.  All  other  parameters  were  fixed  at  literature  values,  or  related  to  other  adjustable 

1 9  parameters  as  depicted  in  Table  2. 

20  2. 2, 4. 3.  Parameter  optimization  and  sensitivity  analysis 

21  PEST  (Parameter  ESTimation),  a  model-independent  parameter  estimator  developed  by 

22  Doherty  (1994).  was  selected  for  parameter  optimization  based  upon  its  extensive  use  and  robust 


16 


]  performance  (e.g.  Doherty.  2003;  Gallagher  and  Doherty.  2006;  Tonkin  and  Doherty,  2008). 

2  PEST  employs  a  Levenberg-Marquardt  nonlinear  least  squares  fitting  algorithm  (Levenberg, 

3  1944;  Marquardt,  1963)  to  adjust  input  parameters  to  minimize  the  objective  function 

4  <j)  =  X(Mz)",  where  w,  and  r,  express  the  weighting  factor  and  the  difference  between  model 

i=i 

5  outcomes  and  the  measurement  for  the  fith  observation  (error),  respectively.  Using  the  initial 

6  values  shown  in  Table  2,  the  calibration  process  searched  for  the  optimal  set  of  biokinetic 

7  parameters  that  resulted  in  the  smallest  difference  between  the  model-predicted  and 

8  experimentally  determined  effluent  concentrations  of  chlorinated  ethenes.  A  total  of  100  column 

9  effluent  concentration  observations  for  PCE,  TCE,  e/s-DCE,  VC,  and  ethene  (20  observations  for 

10  each  of  the  chlorinated  ethenes  and  ethene  distributed  from  5  to  50  days  of  column  operation) 

1 1  were  chosen  for  this  calibration.  To  account  for  the  large  (approximately  1 0  times  greater)  cis- 

12  DCE  concentrations  relative  to  the  other  chlorinated  ethenes,  the  cfy-DCE  concentrations  were 

13  weighted  by  a  0.10  factor.  This  weighting  ensured  that  the  relative  importance  of  each 

14  chlorinated  ethene  was  similar.  Yield  values  were  constrained  to  the  range  reported  in  the 

15  literature  (Table  2).  Fitted  parameters  were  log-transformed  to  approximate  linearity,  which 

16  facilitates  optimization  (Doherty,  1994). 

17  Parameter  optimization  using  PEST  was  distributed  to  multiple  local  network  computers 

18  operated  in  a  parallel  configuration.  The  minimum  number  of  model  runs  for  each  optimization 

19  iteration,  the  number  required  to  calculate  the  Jacobian  matrix  (the  matrix  of  the  derivative  of 

20  observations  with  respect  to  parameters),  was  eight.  The  optimization  process  was  assumed  to 

21  be  complete  when  three  sequential  iterations  failed  to  decrease  the  value  of  the  objective  function. 

22  Sensitivity  of  the  100  observations  to  the  20  biological  model  parameters,  i.e.,  the 

23  maximum  utilization  rates,  yields,  and  half-saturation  constants  (15  parameters),  PCE  and  cis- 
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1  DCE  inhibition  constants  (2  parameters),  and  the  endogenous  decay  rates  of  (he  PER,  GEOSZ, 

2  and  DHC  (3  parameters)  was  examined  using  a  Composite  Sensitivity  (CS)  (Doherty,  1994)  of 

3  each  parameter p,  for  each  specific  group  containing  n  observations  Of. 


5  where  w,  is  the  weighting  factor  for  observation  o,.  The  composite  sensitivity  (CS)  is  a  measure 

6  of  the  contribution  of  the  independent  variable  to  the  variance  of  the  dependent  variable  and 

7  provides  insight  into  the  impact  of  the  input  parameters  on  the  output  of  interest  (observation). 

8  In  this  study,  5  observation  groups:  PCE,  TCE,  ew-DCE.  VC,  and  ethene  effluent  concentrations, 

9  each  consisting  of  20  observations,  were  considered  to  quantify  the  contribution  of  each 

10  parameter  to  the  observed  concentrations.  Additionally,  the  CS  for  all  20  parameters,  based  on  a 

11  single  group  of  all  100  observations,  was  used  to  evaluate  the  overall  importance  of  each 

12  parameter.  The  higher  the  CS  value,  the  more  influence  that  parameter  had  on  the  observations. 

1 3  A  zero  value  indicates  (hat  parameter  variation  had  no  influence  on  the  observation. 

14 

15  3.0  Results  and  Discussion 

16  3. 1  Modeling  Column  Data 

17  Comparisons  of  model  simulations  using  the  optimized  parameter  set  (solid  lines)  with 

18  experimental  measurements  are  shown  in  Figure  2.  For  the  optimal  fit  shown  in  Figure  2  (Table 

19  2  calibrated  parameters),  the  model  was  run  a  total  of  208  times.  This  comparison  reveals  that 

20  the  calibrated  effluent  concentrations  match  the  general  trends  of  the  experimental  measurements 

21  wrell.  Relative  errors  between  experimental  and  simulated  effluent  concentrations  (r,)  were 

22  typically  less  than  30%  for  each  of  the  chlorinated  cthenes  and  ethene  except  at  late  (>45  days) 


18 


1  times  (discussed  below).  As  in  the  experiment,  PCE  dissolution  and  subsequent  dechlorination 

2  occurred  with  minimal  TCE  accumulation.  Accumulation  of  cv.v-DCE  and  VC  was  observed 

3  beginning  near  day  12,  presumably  once  the  microbial  populations  had  become  acclimated. 

4  Comparison  of  simulated  and  experimental  cfv-DCE  concentrations  demonstrated  close 

5  agreement  across  the  observed  concentration  range  (average  r,  <  0.5%)  and  throughout  most  of 

6  the  experimental  period,  However,  comparison  of  simulated  and  observed  VC  and  ethene 

7  effluent  concentrations  revealed  significant  deviations  after  day  45.  Experimental  VC  and 

8  ethene  concentrations  were  over-predicted  immediately  following  the  depletion  of  c/s-DCE 

9  around  day  46.  Given  that  the  predicted  cw-DCE,  VC  and  ethene  concentrations  at  early  times 

10  (<  40  days)  showed  good  agreement  with  the  experimental  data  (average  r,  <  8%),  there  may  be 

1 1  unidentified  mechanisms  (e.g.,  an  unidentified  biotic/abiotic  process  or  loss)  not  captured  by  the 

12  model  that  resulted  in  lower  dissolved  phase  VC  and  ethene  concentrations.  Overall,  in  the 

13  experimental  system,  83%  of  the  initial  PCE  mass  was  recovered  as  cri-DCE,  VC  and  ethene. 

14  Attempts  to  mimic  this  experimental  behavior  in  the  model  by  reducing  the  intial  chloroethene 

1 5  mass  in  the  system  were  unsuccessful,  suggesting  that  some  sort  of  mass  sink  was  present  and 

16  became  more  pronounced  over  time.  One  explanation  of  the  observed  mass  loss  is  the  formation 

17  of  gas-filled  pockets  in  the  experimental  column  (e.g.,  Yang  and  McCarty.  2002;  Amos  et  ah, 

IS  2009).  VC  and  ethene  in  the  plume  region  could  partition  into  these  pockets,  acting  as  a  sink  tor 

19  these  compounds  and  reducing  their  aqueous  phase  concentrations.  Alternatively,  recent  work 

20  has  suggested  that  transformation  product  aqueous  concentrations  may  decrease  due  to  back 

21  partitioning  to  the  NAPL  phase  (Ramsburg  et  ah,  2010a,  2010b).  The  absence  of  kinetic  and 

22  equilibrium  partitioning  data  for  the  chlorinated  daughter  products  and  ethene  with  the  mixed 

23  NAPL,  prevents  quantification  of  this  process.  However,  the  relatively  small  DNAPL  source 
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1  region  (~10  cm)  suggests  any  back  partitioning  is  likely  a  secondary  effect  (estimated  at  <1%  of 

2  the  total  chloroethene  mass  assuming  equilibrium  partitioning  data  from  Ramsburg  et  al.,  2010a). 

3  While  the  model  described  effluent  concentration  data  quite  well,  a  more  revealing 

4  indicator  of  model  performance  is  the  ability  of  the  calibrated  model  to  predict  the  microbial 

5  mass  and  transformation  product  measurements  in  the  source  zone  and  the  plume  region,  data 

6  that  were  not  used  in  the  calibration  process.  A  comparison  of  predicted  and  measured  data  from 

7  side  port  samples  along  the  length  of  the  column  provides  an  independent  assessment  of  the 

8  model’s  capability  to  simulate  the  reductive  dechlorination  process. 

9  PCE,  c/s-DCE.  and  VC  concentrations  collected  from  side  ports  along  the  length  of  the 

10  column  on  days  16,  23  and  42  were  compared  to  corresponding  simulated  concentrations  in 

1 1  Figure  4.  At  day  1 6  (Figure  4a),  the  calibrated  model  produced  reasonable  predictions  of  cis- 

12  DCE  and  VC  concentrations  along  the  length  of  the  column  (r,  typically  less  than  10%  for  cis- 

13  DCE  and  20%  for  VC);  however,  PCE  concentrations  were  significantly  under- predicted  (r,  ~ 

14  70%).  This  result  suggests  that  simulated  microbial  activity  in  the  source  zone  was  over 

15  predicted  at  day  16  or  that  the  Sherwood  number  employed  in  the  model  under-estimated  the 

16  mass  transfer  rate,  leading  to  the  under-estimation  of  dissolved-phase  PCE  concentrations  in  the 

17  source  region  (x  <  0.15m).  No  attempt  wras  made  to  adjust  the  mass  transfer  rate  since  the 

18  dissolution  kinetic  correlations  employed  in  this  work  are  well  established  (Powers  et  ah,  1994). 

19  In  the  day  23  comparison  (Figure  4b),  PCE  was  still  present  in  the  source  zone  and  significant 

20  amounts  of  cA-DCl:  and  VC  had  formed  in  the  experimental  system.  Model  simulations 

2 1  provided  accurate  predictions  of  experimental  observations  of  cri-DCE  (r,  typically  less  than 

22  1 0%)  and  predictions  within  an  order-of-magnitude  for  VC  (r,  typically  less  than  70%),  but  again 

23  tended  to  under-estimate  the  PCE  concentration  in  the  source  region.  When  compared  to  day  42 
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1  experimental  results  {Figure  4c).  the  model  captured  the  general  behavior  of  the  cu-DCE  and  VC 

2  concentrations  along  the  length  of  the  column  (average  rcis.DCE  <  5%  and  rvc  typically  less  than 

3  60%).  The  variability  in  the  cA-DCE  experimental  measurements  along  the  length  of  the  column 

4  is  likely  due  to  the  spatial  variability  in  PCE  saturation  as  the  PCE  was  depleted  from  the  source 

5  zone  at  later  times.  Given  that  the  model  is  one-dimensional  (1-D),  it  is  not  capable  of  capturing 

6  small  scale  nonuniformity  in  domain  conditions  that  emerge  over  time,  although  it  does  capture 

7  the  general  increase  and  plateau  of  cis- DCE  and  VC  concentrations. 

8  In  addition  to  predicted  transformation  product  concentrations,  model  predictions  of 

9  microbial  growth  were  evaluated,  a  comparison  that  is  generally  absent  in  the  literature. 

10  Unattached  (“planktonic”)  aqueous  phase  microbial  cells  were  measured  from  the  side  ports 

11  during  column  operation  (Figure  5).  The  model  assumes  a  total  active  biomass  consisting  of 

12  both  attached  and  unattached  cells.  To  facilitate  comparisons  between  experimental  biomass 

13  measurements  (unattached)  and  the  model  predictions  (total),  the  ratio  between  total  and 

14  unattached  biomass  was  estimated  at  16,  23,  and  42  days  (Table  3).  This  ratio  was  assumed  to 

1 5  be  a  constant  in  the  source  zone  (x  <  0. 1 5  m)  and  the  plume  region  (0.15  m  <  x  <  0.60  m)  at  each 

16  sampling  time,  and  was  adjusted  within  the  experimental  range  given  by  the  initial  (day  0)  and 

17  final  (day  42)  measurements  to  give  the  best  visual  lit  to  the  GEO  and  DIIC  spatial  distributions 

18  along  the  column.  Figure  5  depicts  experimental  and  simulated  unattached  GEO  and  DIIC 

1 9  biomass  data.  The  simulated  biomass  profiles  scaled  to  reflect  only  the  unattached  biomass  are 

20  in  agreement  with  experimental  measurements  at  most  of  the  side-ports  on  day  16  (median  rGEO 

21  =  10.5%  and  i'dhc  =  7.0%),  Although  the  model  over-estimated  the  DHC  biomass  in  the  plume 

22  region  by  day  23  and  day  42,  model  predictions  are  still  within  approximately  one  order-of- 

23  magnitude  of  experimental  observations  (Figure  5b,  c;  median  rGE0  =  20%  and  roue  =  170%  at 
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1  day  23  and  rGEO  ~  35%  and  rDHC  =  112%  at  day  42).  Note  that  fitting  the  total  to  unattached 

2  biomass  ratio  served  to  scale  (i.e.,  to  adjust  the  simulated  curve  up  and  down)  the  spatial 

3  distribution,  but  did  not  alter  the  spatial  trends  in  biomass  along  the  column  length.  Thus,  the 

4  mathematical  model  served  as  a  good  predictor  for  the  observed  microbial  growth  patterns, 

5  Although  adjustment  of  kinetic  parameters  (e.g.,  maximum  growth  rates  or  yield  coefficients) 

6  may  have  provided  an  equally  reasonable  fit  to  the  effluent  concentrations,  simulations  based 

7  upon  modified  microbial  kinetic  parameters  failed  to  adequately  reproduce  the  side  port  data, 

8  suggesting  that  the  fitted  parameters  were  unique  and  provided  the  best  prediction  of  effluent  and 

9  side-port  experimental  results. 

10 

1 !  3.2.  Sensitivity  Analysis 

12  3.2.  /.  Composite  sensitivity  of  microbial  kinetics 

13  Based  upon  the  results  presented  in  the  previous  section,  the  calibrated  model  was  used  to 

14  assess  the  sensitivity  of  reductive  dechlorination  activity  to  model  parameters  by  identifying  the 

1 5  relative  importance  of  each  of  the  kinetic  parameters  (Table  2)  using  the  CS.  Figure  6  depicts 

16  the  CS  for  each  parameter  for  each  group  of  concentration  measurements.  Figures  6a  and  6b 

17  indicate  the  first  two  reductive  dechlorination  steps  from  PCE  via  TCE  to  cE-DCE  are  most 

1 8  sensitive  to  the  PCE  maximum  utilization  rate  ( kff  ),  i.e.,  that  the  accumulation  of  cis- DCE  was 

19  influenced  by  the  degradation  of  the  higher  chlorinated  contaminants  more  so  than  the  rate  of 

20  c/.v-DCE  dechlorination  (low  CS  lor  kmaxDCE).  The  cA-DCE-to-VC  dechlorination  step  appears  to 

21  be  nearly  as  sensitive  to  the  growth  of  the  GEO  population,  as  indicated  by  the  high  CS  for  Ypqe 

22  (Figure  6c).  The  degradation  of  cA-DCE  ( kff ),  as  well  as  the  growth  of  the  DFIC  population 
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1 


( Y!H  /; ,  kjH<  ),  controlled  the  second  phase  of  the  dechlorination  process  (i.e.,  c/s-DCE  to  ethene, 

2  Figure  6d  and  6e).  Interestingly,  Figure  6f  shows  that  the  overall  model  results  are  most 

3  sensitive  to  the  m-DCE  maximum  utilization  rate  ( k^J  )  and  yield  of  cis- DCE  (  Y!K  h ),  followed 

4  by  the  parameters  governing  PCE  degradation.  This  finding  suggests  that  PCE  and  c/s-DCE 

5  degradation  are  the  key  steps  in  the  sequential  dechlorination  process.  Apparently,  once 

6  dechlorination  from  c/s-DCE  to  VC  occurs,  the  final  dechlorination  step  to  ethene  is  not  limiting 

7  in  this  system.  Sensitivities  of  kinetic  parameters  associated  with  lactate  fermentation  (k^, 

8  K's"‘ ,  )  were  not  significant  in  this  modeling  application  due  to  the  consistent  delivery  of 

9  sufficient  electron  donor  to  fully  support  reductive  dechlorination  reactions.  These  results  are 

10  consistent  with  earlier  findings  of  Kouznctsova  et  al.  (2010),  who  showed  that  dechlorination  in 

1 1  a  batch  system  was  most  sensitive  to  ( )  and  ( Y,x  A ).  However,  their  system  did  not  consider 

1 2  the  PCE-to-TCE  dechlorination  step  and  did  not  include  NAPE  dissolution. 

13  3.2.2.  Inhibition 

14  Although  the  above  sensitivity  analysis  indicated  that  model  predictions  were  not 

15  sensitive  to  the  cis- DCE  inhibition  coefficient  (  ks*h:  )  for  these  experimental  conditions. 

1 6  inhibitory  effects  of  polychlorinated  ethenes,  and/or  toxicity  due  to  the  decrease  in  pH  during  the 

17  dechlorination  process  are  often  cited  as  a  reason  for  incomplete  dechlorination  to  ethene  (e.g., 

18  Adamson  et  ah,  2004:  Haest  et  ah,  2010;  Sabalowsky  and  Semprini,  2010a,  2010b;  Amos  ct  ah, 

19  2008,  2009).  The  calibrated  model  provides  a  platform  for  a  more  detailed  analysis  to  assess  the 

20  influence  of  inhibition  on  the  complete  dechlorination  of  cis- DCE  and  VC  to  ethene.  Figure  7 

21  depicts  the  simulated  accumulation  of  c/.v-DCE,  VC,  and  ethene  in  the  column  effluent  using  a 

22  no-inhibition  scenario  and  the  inhibition  coefficients  listed  in  Table  2.  In  the  former  case,  total 
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1  ethene  production  increased  ten-fold  and  total  m-DCE  and  VC  formation  were  reduced  by  25% 

2  and  20%,  respectively.  Thus,  while  the  model  results  were  relatively  insensitive  to  the 

3  magnitude  of  the  inhibition  coefficient,  neglecting  inhibition  could  lead  to  substantial 

4  underestimation  of  the  accumulation  of  VC  and  overestimation  of  the  complete  conversion  to 

5  ethene, 

6  3.2.3,  Column  residence  time 

7  The  experimental  results  of  Amos  et  al.  (2009)  suggest  that  the  column  residence  time  {ca. 

8  2.8  days)  was  insufficient  for  complete  dechlorination,  and  thus,  a  series  of  numerical  studies 

9  was  performed  to  examine  the  influence  of  residence  lime  {i.e,,  column  length  and  aqueous 

10  phase  (low  rate)  on  simulation  results.  In  these  numerical  studies,  the  column  length  was  first 

1 1  halved,  and  then  doubled,  relative  to  the  initial  column  length  (60  cm),  while  all  other  parameters 

12  were  kept  constant  (see  Tables  1  and  2).  Simulation  results  for  the  shortened  column  resulted  in 

13  a  decrease  in  the  total  production  of  cis- DCE,  VC,  and  ethene  by  19%,  53%,  and  73%, 

14  respectively  (Figure  8).  As  expected,  the  reduced  column  length  resulted  in  a  residence  lime  (ca. 

15  1.4  days)  that  was  insufficient  for  complete  PCE  conversion  to  c/.v-DCE.  When  the  column 

16  length  was  doubled,  less  cw-DCE  accumulated,  while  VC  and  ethene  production  increased 

1 7  significantly.  The  42%  reduction  in  c/s-DCE  production  (up  to  day  55)  in  the  longer  column 

1  8  indicated  that  the  original  column  length  was  insufficient  to  achieve  complete  dechlorination  to 

19  ethene.  In  the  longer  column  with  increased  residence  time  (ca.  5.6  days),  the  DHC  was 

20  predicted  to  dechlorinate  efv-DCE  more  effectively,  resulting  in  a  4-fold  and  a  9-fold  increase  in 

21  VC  and  ethene  production,  respectively  (Figure  8). 

22  A  simpler  means  to  modify  residence  time  in  the  column  is  to  change  the  aqueous  phase 

23  flow  rate.  However,  changing  the  aqueous  phase  flow  rate  not  only  affects  column  residence 
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time,  but  also  impacts  the  rate  of  PCE  mass  transfer  from  the  NAPL  to  the  aqueous  phase  (see 
Eq.5).  In  this  set  of  simulations,  the  aqueous  phase  flow  rate  was  first  halved  and  then  doubled 
for  comparison  with  the  original  simulation  conditions.  As  expected,  maximum  PCE  discharge 
(pg/min)  during  the  initial  4  days  was  linearly  proportional  to  the  flow  rate  due  to  the  fact  that 
the  equilibrium  aqueous  phase  PCE  concentration  was  reached  in  all  three  cases  (Figure  9). 
When  the  flow  rate  was  doubled  (residence  time  of  ca.  1.4  days),  c/i-DCE  and  VC  reached 
higher  maximum  discharge  levels  after  the  flow  interruption  (dashed  vertical  lines  in  Figure  9), 
but  were  also  depleted  faster  (Figure  9)  than  in  the  base  case  simulation,  resulting  in  less  total 
mass  recovered  (Figure  10).  Ethenc  production  occurred  earlier  in  the  doubled  flow  rate 
simulation,  but  decreased  earlier  due  to  less  VC  available.  In  the  reduced  flow  rate  simulation 
(residence  time  of  ca.  5.6  days),  VC  degradation  was  inhibited  by  the  elevated  cw-DCE 
concentrations  (see  eq.  7.  k/KE  =  0.52.4  //g/L,  Table  2)  resulting  in  minimal  cthene  production 
by  the  end  of  the  simulation  (Figure  10).  Thus,  despite  the  longer  residence  time,  microbial 
inhibition  resulting  from  elevated  levels  of  cis-DCE  prevented  complete  dechlorination  to  cthene. 

Interestingly,  the  model  predictions  indicate  that  the  total  mass  of  chlorinated  ethenes  and 
ethene  recovered  in  the  column  effluent  increases  in  proportion  to  the  flow  rate  due  to  the 
increase  in  PCE  dissolution  (Figure  1 1 ).  At  the  high  flow  rate,  more  PCE  dissolved  into  the 
aqueous  phase  and  was  recovered  in  the  effluent,  leading  to  a  higher  total  mass  recovery. 
However,  less  PCE  was  dechlorinated  to  c/.v-DCE  (45%  less  than  baseline),  VC  (39%  less),  and 
ethene  (73%  less).  Although  91%  and  169%  more  c/.v-DCE  and  VC,  respectively,  were 
produced  under  the  slow  flow'  condition,  PCE  dissolution  was  reduced,  resulting  in  greater  PCFI 
mass  remaining  in  the  NAPE.  Additionally,  38%  less  ethene  was  produced  in  the  low  flow 
column  due  to  inhibition  of  VC  reductive  dechlorination  due  to  the  high  c/.v-DCE  concentrations. 
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1  These  findings  indicate  that  residence  time  in  the  bioactive  zone  is  a  key  factor  determining 

2  detoxification  (i.e.,  cthene  formation).  The  refined  model  will  be  useful  to  optimize  the  design 

3  of  source  zone  bioremediation  in  order  to  achieve  treatment  objectives,  whether  that  means 

4  removal  of  contaminant  mass  from  the  system  as  rapidly  as  possible,  or  complete  dechlorination 

5  to  non-toxic  ethene,  or  both. 

6 

7  4. 0  Conclusions 

8  A  multiphase  compositional  simulator  was  designed  to  accurately  represent  the  sequential 

9  reductive  dechlorination  process  converting  PC'E  to  ethene  by  a  microbial  consortium  containing 

10  Geobacter  and  Dehalococcoides  populations  under  dynamic  conditions.  The  dechlorination 

1 1  process  resulted  in  enhanced  dissolution  of  PC  Li  from  a  source  zone  containing  a  mixed  N  A  PL  at 

12  residual  saturation.  The  model  was  calibrated  to  effluent  chlorinated  ethenes  and  ethene 

13  concentrations  measured  during  a  bioenhanced  dissolution  column  experiment  by  optimizing  the 

14  eight  most  sensitive  microbial  kinetic  parameters  using  a  best-fit  calibration  process.  Model 

15  simulations  generally  captured  measured  effluent  concentrations  of  all  five  chlorinated  ethene 

16  constituents  (predicted  to  experimental  relative  errors  less  than  30%).  Simulated  VC  and  ethene 

17  concentrations  were  higher  than  observed  values,  possibly  due  to  gas  evolution  and  gas  phase 

18  partitioning,  which  were  not  considered  in  the  model.  Following  calibration,  model  predictions 

19  were  compared  to  concentrations  of  chlorinated  ethenes  and  biomass  measured  in  side-port 

20  samples  collected  at  three  different  times.  This  novel  validation  methodology  demonstrated 

21  model  accuracy  and  applicability  using  a  single  set  of  experimental  results.  This  level  of  model 

22  validation  has  not  been  conducted  for  any  dechlorinating  consortium  lo  date,  and  thus,  the 
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1  present  work  provides  a  unique  opportunity  to  verify  the  conceptual  models  of  microbial  kinetics 

2  commonly  employed  in  dechlorination  simulators. 

3  Using  the  calibrated  model,  a  detailed  sensitivity  analysis  was  performed  to  investigate 

4  the  relative  influence  of  20  microbial  kinetic  parameters  on  the  overall  model  performance.  The 

5  composite  sensitivity  assessment  suggested  that  degradation  of  PCE  and  c/'.v-DCE  are  key 

6  processes  controlling  complete  dechlorination  to  ethene.  A  step-by-step  evaluation  of 

7  dechlorination  kinetics  revealed  that  the  higher  chlorinated  contaminants  (PCE  and  TCE)  were 

8  most  sensitive  to  the  PCE  maximum  utilization  rate,  while  the  lesser  chlorinated  contaminants 

9  (c/a-DCE  and  VC)  were  most  sensitive  to  c/j-DCE  dechlorination  and  DHC  growth  kinetics. 

10  Simulations  that  neglected  inhibition  resulted  in  a  10-fold  increase  in  ethene  production, 

1 1  demonstrating  the  importance  of  inhibitory  factors  on  bioactivity.  Altering  the  column  length  or 

12  flow  rate  to  increase  or  decrease  residence  time  also  indicated  that  complete  dechlorination  to 

13  ethene  was  influenced  by  a  balance  between  sufficient  residence  time  to  achieve  complete 

14  dechlorination  without  allowing  for  accumulation  of  c/5- DC E,  which  inhibits  the  final 

15  dechlorination  step.  Assuming  adequate  parameterization  of  the  bioactive  zone,  the  verified 

16  mode!  may  be  used  to  investigate  bioenhanced  dissolution  in  field  settings,  and  as  a  guide  to 

17  design  biostimulation  and  bioaugmentation  strategies  for  source  zone  treatment  in  conjunction 

1 8  with.  or  following,  aggressive  mass  removal. 
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Figure  captions: 

Figure  1.  Schematic  diagram  of  the  experimental  column  system  (adapted  from  Amos  et  aL, 
2008) 

Figure  2.  Comparison  between  simulated  (solid  line)  and  experimental  (dashed  line  with  marker) 
organic  effluent  concentrations.  The  flow  rate  was  interrupted  between  two  vertical  dashed  lines, 
and  was  decreased  from  0.25  mL/min  to  0.1  mL/min  after  the  vertical  dash-dot  line  on  day  24. 

Figure  3.  Side-port  samples  of  initial  Geobacter  (GEO)  and  Dehalococcoides  (DHC) 
concentrations  along  the  column.  The  initial  DHC  values  at  ports  1  and  2  (5  and  10  cm  -  first 
two  data  points)  are  replaced  by  those  at  day  1 6  (open  symbols). 

Figure  4.  Comparison  of  predicted  (solid  line)  and  experimental  (symbols)  PCE,  f  /.v-DCE,  and 
VC  concentrations  obtained  from  the  side-port  samples  along  the  length  of  the  column  at  (a)  16 
days,  (b)  23  days,  and  (c)  42  days. 

Figure  5.  Comparison  of  predicted  (solid  line)  and  experimental  (symbols)  planktonic  biomass 
(' Geobacter  (GEO)  and  Dehalococcoides  (DHC))  obtained  from  side-port  samples  along  the 
length  of  the  column  at  (a)  16  days,  (b)  23  days,  and  (c)  42  days. 
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Figure  6.  Composite  sensitivity  (CS)  of  kinetic  parameters  (see  Table  2  for  kinetic  parameters) 
to  (a)  PCE,  (b)  TCE,  (c)  cfs-DCE,  (d)  VC,  (e)  ethene,  and  (f)  all  observations. 

Figure  7.  Comparison  of  simulated  cumulative  production  of  c/s-DCE,  VC,  and  ethene  for  the 
base  case  simulation  (open  symbols)  and  no-inhibition  simulation  (closed  symbols). 

Figure  8.  Simulated  total  mass  (kg)  of  cix- DCE,  VC.  and  ethene  produced  by  day  55  for  the  base 
case,  and  for  cases  with  the  column  length  (L)  reduced  by  a  factor  of  0.5  or  increased  by  a  factor 
of  2. 

Figure  9.  Mass  effluent  flux  (pg/min)  of  PCE,  c/s-DCE,  VC,  and  ethene  for  the  base  case  (solid 
line,  no  symbols),  and  for  cases  with  the  Darcy  velocity  (qw)  reduced  by  a  factor  of  0.5  (dashed 
tine  with  open  symbols)  or  increased  by  a  factor  of  2(dashed  line  with  closed  symbols). 

Figure  10.  Simulated  total  mass  (kg)  of  cfs-DCE,  VC,  and  ethene  produced  by  day  55  for  the 
base  case,  and  for  cases  with  the  Darcy  velocity  (qv)  reduced  by  a  factor  of  0.5  or  increased  by  a 
factor  of  2. 

Figure  11.  Simulated  total  chloroelhene  and  ethene  mass  left  in  the  column  (left-hand  y-axis, 
open  symbols)  and  recovered  (right  hand  y-axis,  closed  symbols)  for  the  base  case  (diamonds), 
and  for  cases  with  the  Darcy  velocity  (qw)  reduced  by  a  factor  of  0.5  (circles)  or  increased  by  a 
factor  of  2  (triangles). 
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Table  1.  Physical  Column  Experimental  Parameters 


Parameters 

Value 

Porosity,  - 

0.374a 

Median  grain  size  d;o,  mm 

0.32a 

Uniformity  index  (d6o/  dio),  - 

1.86b 

Bulk  soil  density,  g/cm3 

2.03° 

Intrinsic  permeability,  darcy 

42. 0a 

Dispersivity,  m 

0.01d 

PCE-Hexadecane  (HD)  mole  ratio 

0.25/0.75 

PCE/HD-NAPL  density,  g/cm3 

0.863 

Bottom  water  flux,  mL/min 

0.25.0.0.1s 

Column  length,  cm 

60 

Column  diameter,  cm 

4.8 

Total  pore  volume,  cm3 

407c 

Source  zone  position  &  saturation 

0-10  cmr& 0.21s 

Transition  zone  position  & 
saturation 

10-15  cm r&  0.21  to  0.08g 

3  Estimated  from  Suchomel  et  aL  2007, 
b  Estimated  from  Santos  and  Barros.  2010. 
c  Estimated  from  Amos  el  al.,  2009. 
d  Assumed. 

c  Initial  flow  0,25mL/min;  flow  interruption  at  19th  day  for  22  hours;  0.1  mL/min  after  24th  day, 
f  Visual  observation 

g  Total  NAPL  volume  was  14  ml  based  on  Amos  et  ahf  2009. 
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Table  2,  Initial  and  calibrated  model  coefficients  for  dechlorination  used  in  simulations 


Parameters 

Type  Literature3 

Initial2 

Calibrated 

Maximum  utilization  rate  klmAX  (ptmol/mg-day) 

Lactate 

Fixed  206 

206 

- 

PCE 

Adjustable  1.4-117 

14 

94 

TCE 

Adjustable  2.4  -  366 

24 

80 

DCE 

Adjustable  1.7-48 

18 

30 

VC 

Adjustable  2.6-48 

3.2 

6.9 

Yield  YL  (g  cell/mol  substrate  )  ~ 

Lactate 

Fixed  1. 5-6.3 

3.9 

- 

PCE 

Adjustable  7.54-22,6 

15.1 

17.9 

TCE 

Tied 

15.1 

17.9 

DCE 

Adjustable  5.89-8.25 

7.07 

5.89 

VC 

Tied  =Y,k'k 

7.07 

5.89 

Half  saturation  constant  K"s  (pM) 

Lactate 

Fixed  2.5 

2.5 

- 

PCE 

Fixed  0.11-2.8 

0.54 

- 

TCE 

Fixed  0,54-1.5 

0.54 

- 

DCE 

Fixed  0,54-3.3 

0,54 

- 

VC 

Fixed  0.54  -  360 

0.54 

- 

Inhibition  coefficients  k\  (gM/ 

PCE 

Fixed  =  K!’ce 

0.54 

- 

DCE 

Fixed  =  k'sc 

0.54 

- 

Microbe  decay  rate  kd  (1  /day) 

Fermentor 

Fixed  0.050 

0.050 

- 

Geobacter 

Adjustable  0.050 

0.050 

0.063 

Dehaio  cocco  ides 

Adjustable  0.050 

0.050 

0.048 

1  Literature  values  for  the  Maximum  utilization  rate  and  half  saturation  constant  were  taken  from 
Clapp  et  al,  (2004)  and  Lee  et  ak  (2004);  Yield  for  lactate  and  chlorinated  ethenes  are  taken  from 
Sedtger  et  at. (2002)  and  Amos  et  al.  (2009). 

2  Initial  values  for  calibration  process  are  estimated  based  on  the  literature. 

and  ^  are  tied  to  *  n  and  ^ tKL  respectively,  and  are  not  included  in  the  optimization 
process. 

4  Inhibition  coefficients  are  tied  to  the  corresponding  half-saturation  constant  as  discussed  in 
Christ  and  Abriola  (2006). 
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Table  3:  Total/unattached  ratio  of  biomass  in  column 


Source  zone 

Plume  zone 

Initial1' 

Geo  hue  ter 

20 

20 

Dehalococcoides 

20 

20 

Day  16b 

Ge  abac  ter 

20 

10 

Dehalococcoides 

20 

10 

Day  23" 

Geobacter 

20 

10 

Dehalococcoides 

10 

5 

Day  42h 

Geobacter 

20 

2 

Dehalococcoides 

10 

2 

Final" 

Geobacter 

4-20 

1 

Dehalococcoides 

10 

1.25 

3  Ratio  based  on  experimental  measurements  prior  to  column  operation  and  at  completion 
of  the  experiment. 

h  Ratio  based  on  best-fit  to  side  port  samples  of  biomass  concentrations. 
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